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ABSTRACT 



Radio observations of galaxy clusters show that there are fxG magnetic fields 
permeating the intra-cluster medium (ICM), but it is hard to accurately constrain 
the strength and structure of the magnetic fields without the help of advanced 
computer simulations. We present qualitative comparisons of synthetic VLA ob- 
servations of simulated galaxy clusters to radio observations of Faraday Rotation 
Measure (RM) and radio halos. The cluster formation is modeled using adap- 
tive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with 
the assumption that the initial magnetic fields are injected into the ICM by ac- 
tive galactic nuclei (A GNs) at high redshift. In addition to simulated clusters in 



Xu et al.l (120101 . 12OIII ). we present a new simulation with magnetic field injections 



from multiple AGNs. We find that the cluster with multiple injection sources is 
magnetized to a similar level as in previous simulations with a single AGN. The 
RM profiles from simulated clusters, both \RM\ and the dispersion of RM {(Tjim), 
are consistent at a first-order with the radial distribution from observations. The 
correlations between the a^M and X-ray surface brightness from simulations are 
in a broad agreement with the observations, although there is an indication that 
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the simulated clusters could be slightly over-dense and less magnetized with re- 
spect to those in the observed sample. In addition, the simulated radio halos 
agree with the observed correlations between the radio power versus the cluster 
X-ray luminosity and between the radio power versus the radio halo size. These 
studies show that the cluster wide magnetic fields that origi nate from AGNs and 
are then amplified by the ICM turbulence (jXu et al.l 120101 ) match observations 
of magnetic fields in galaxy clusters. 



Subject headings: galaxies: clusters: general — large-scale structure of Universe- 
magnetic fields — methods: numerical — MHD 



1. Introduction 



The detections of lar ge-scale radio synchro t ron emissions from galaxy clusters, called 
radio halos and relics (see 



Carilli & Tavlor 


2002; 


Ferrari et al. 


2008; 


Giovannini et al. 


2009 



Feretti et al.l 120121 ). indicated that the intra-cluster medium (ICM) of galaxy clusters are 
permeated with cluster wide magnetic fields. The radio halos are diffusively extended over 
~ 1 Mpc, covering the whole clusters, while the radio relics are usually at the edges of the 



clust e rs as long filameiits and could be related to recent cluster mergers (e.g. iBagchi et al. 
20061 : Ivan Weeren et al.l |2010| ) . By assuming that the magnetic energy is comparable to the 
total energy in relativistic electrons, one often deduces that the volume-averaged magnetic 
fields in the cluster halos are ~ 0.1 to 1.0 /iG and the total magnetic energy can be as high 



as 10^^ erg (IFerettil Il999l ) in the ICM of a galaxy cluster. Recently, observations of radio 
halos have been used to study the structure of the cluster wide magnetic fields b y comparing 



obser vations with mock halos from turbulent magnetic fields by construction (IVacca et al. 
2010h . 



The Faraday rotation measure (RM) is another importa nt observational techn ique to 



study the pro perties of magnetic fi e lds in galaxy clus t ers (s e e ICarilli fc Tav^ 



oi 



20021). Stud- 



ies of RM by lEilek fc OwenI fl2002f ): iTavlor fc PerlevI fll993h : ICokate fc Lil ((20001) have sug- 
gested that the coherence scales of magnetic fields can range from a few kpc to a few hun- 
dred kpc. Recently, the patchy RM maps were further studied to predict that the ICM 
magnetic fields have a Kolmogorov-like turb ulent spectrum in the cores of some clusters 
(lEnfilin fc Vogtll2006l : iKuchar fc EnBlinll201ll ). RM is wi dely used to estimate the strengths 



of th e ICM magnetic fie l ds from the cluster centers (e.g. iTavlor fc Perleyl Il993l : iTaylor et al. 



2001; Laine et al. 


2008 


ters (e.g.. 


Clar 


ffi et al. 


Bonafede et al. 


2010 





Guidetti et al 



2001 



2010 



Murgia et al 



Vacca et al.l 



2004 



2012!) to the outer parts of clus- 



Govoni et al. 2006; Guidetti et al. 2008 



20101). Though there are large uncertainties in measuring the magnetic fields. 
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these studies all suggest that the strengths of magnetic fields in the cluster centers range 
from 10s yuG in the cool core clusters to a few fiG in the non-cool core clusters, and the 
magnetic fie l d stre ngths drop gradually to sub /iG at the edges of clusters. In addition, 
Dolag et al.l ( 120011 ) showed that there is a power law correlation between two observables, 
the dispersions of RM {aRu) and the X-ray surface brightness (Sx), and suggested that the 
strengths of magnetic fields are related to the gas tempera t ures o f the ICM. This aRM-^x 
correlation was further confirmed recently by iGovoni et al.l (120101 ) with more observational 
data. 

Though the ICM magnetic fie lds may not be dynamically important for the cluster 
formation as a whole (IWidrowll2002l ). the presence of magnetic fields and their structure are 
important to many interesting astrophysics and plasma physics phenomena, e. g., anisotropic 



thermal conduction and acceleration and propaga tion of cosmic rays (e.g., iParrish et al. 
20091 : iBrunetti fc LazarianI 120071 : lEnBlin et al.l |20 111 ). Some of these processes, like thermal 
conduction, may p lay an imp ortant role in the dynamics and thermal equilibrium in cores 
of galaxy clusters (IVoitll201ll ). 



Since the formation of galaxy clusters and evolution of the ICM magnetic fields are 
highly non-linear, thus are too complicated to be studied analytically, numerical simulations 
have been playing an i r nportant part iri studying the thermal properties of the ICM (e.g., 

'201l[). and the ICM 



Brvan fc NormanlllQQsl : iMotl et ahlbooi iNaeai et al.lboOTi: IVazza et al 



magnetic field evolut ion (e.g., bolag et all l2002[ iDubois fc Teyssierl 



2008 



Xu et al. 2009 



Bonafede et al.ll201lf ) of galaxy clusters. Although the existence of cluster-wide magnetic 
fields is well accepted, their origin, which may be import ant to the evolution of magneti c 
fields during the course of cluster formation, is still unclear ( 1Widrowll2002l : iDolag et al.ll2008l ) . 
MHD cluster formation simulations have been performe d with different initial magnetic fields, 
including random or unif orm fields from high redshifts ( Dolag et al. 2002:lDubois fc Tevssier 
20081 : IPubois et alilgoggi) or fr om the outfiows of normal galaxies (IDonnert et al.ll2009l ) or 
active galaxies ( Xu et alTboOol ). The cluster magnetic fields of all these simulations at low 
redshifts are roughly in agreement with each other, with /xG in the cluster centers and 
decreasing with radius, and can produce some mock observations similar to available radio 
observations. So the exact origins of magnetic fields may not be differentiated by the cluster 
magnetic fields at the current epoch. 

Large scale magnetized radio jets and lobes from AGNs serve as a very intriguing source 
of seed magnetic fields in clusters, because they could carry l arge amounts of magnetic energy 



and fiuxes and distribute them to scales of hundreds kpc (iBurbidgd Il959l : iKronberg et al. 



200ll : ICroston et al.l l2005l : iMcNamara &: NulsenI 120071 ). The simulation presented in Xu et 
al (2009) suggests that cluster-wide [iG magnetic fields can be produced from a single AGN 
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that is further am plified by the sraall s cale dynamo ( iBrandenburg Sz Subramaniaru l2005l ). 
Additional studies (IXu et al.ll2010U201ll ) show that this mechanism is mostly insensitive to 
the injected energy and redshift, as well as cluster mass and formation history. The RM 
from these simulations show similar features of long filaments and small scale bands as in 
the observations from VLA. But to know how the AGN magnetic fields fit the real cluster 
magnetic fields, detailed comparisons between the synthetic observations of simulated clus- 
ters and the radio observations are needed. Such studies are necessary before applying these 
simulated magnetic fields to make predictions for the next generation of radio observations, 
such as EVLA, LOFAR, and SKA, or to the study of the particle acceleration by MHD 
turbulence in the ICM. 

In this paper, we first present the results of a new simulation with magnetic field in- 
jections from multiple AGNs. In this simulation, the initial magnetic fields are from 30 
AGNs in a volume that eventually forms a single cluster. This setup of injections is more 
realistic than previously simulations with only one AGN. We then compare the VLA fil- 
tered synthetic radio observations, including RM and radio halos from our high resolution 
MHD cosmological simulations to radio observations. To have a large sample of simulated 
clust e rs, th e comparisons include the simulated clusters from our previous studies (IXu et al. 



2010 



20111). We co mpare the radial profiles of ab solute value of RM ( \RM\) with results in 
Clarke et al.l (1200 ll ). t he (Trm with o b serva tions in lGovoni et al.l ( 120101 ) . and the aRM vs. Sx 



relations with those in iGovoni et al.l (|2010|). For radio halos, we follow the method stated in 
Murgia et al. (j2004 ) and Murgia et al. ( 2009 ) to generate synthetic radio halos by assuming 
energy equipartition between magnetic fields and non-thermal electrons. We then compare 
the statistical properties of the mock radio halos with the real ones. 

The data processing and analysis of the simulations are performed using vl ^ ( Turk et al. 



20M). 



2011) with additional modules for magnetic fields, and the FARADAY code (IMurgia et al. 



This paper is structured as follows. We begin with a description of our new simulation 
in Section [2j In Section [3] , we present the comparisons of synthetic observations with radio 
observations, including RM and radio halos. Finally, we briefly discuss our results and 
conclusions in the last section. We refer to the Appendix for the details of the procedures we 
follow to generate the synthetic observations of RM and radio halos from the simulations. 



^http://yt. enzotools.org 
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2. Simulation of Galaxy Cluster Formation with Multi-AGN Magnetic Field 

Injections 

2.1. Numerical Method and Model 



We first present our new self-consistent high-resolution cosmological MHD cluster for- 
mation simulation with magnetic fields initially from multiple AGNs. The simulation is per- 
formed using cosmo logical MHD code with adaptive mesh refinement (AM R) ENZO-fMHD 



Collins et al.l I20101 ) . Th is code uses the A MR algor i thms developed by iBerger fc Colella 



1989h and lSalsaraj (120011). the MHD solve r of iLi et al.l ( 120081 ). and the constrained transport 



(CT) method of ICardiner fc Stond (120051 ). The simulation here uses an adiabatic equation 
of state, with the ratio of specific heat, T = 5/3, and does not include heating and cooling 
physics or chemical reactions, which are not important in this study. 

The initial c onditi ons of the simulation are generated at redshift z = 30 from an 
Eisenstein fc Hul (119991 ) power spectrum of density fluctuation in a ACDM universe with 
parameters h = 0.73, flm = 0.27, Qb = 0.044, I^a = .73, and = 0.77. These parameters 
are close to the values from WMAP3 observations ( ISpergel et al.l 120071 ). The simulation 
volume is (128 h^^Mpc)^, and it uses a 128^ root grid and 2 level nested static grids in 
the Lagrangian region where a cluster forms. This gives an effective root grid resolution of 
512^ cells (~ 0.35 Mpc) and dark matter particle mass resolution of 1.34 x lO^Mf:^, which 



is 8 t imes better than the dark matter resolution in previous simulations (jXu et al.l 12010 



20111). In this paper, we only study the magnetic field injection and evolution in the central 
cluster out of lots of galaxy clusters formed in the simulation domain. During the course of 
simulation, 7 levels of refinements are allowed beyond the root grid, for a max i mum spatial 
resolution of 7.8125 h^^ kpc, the same as the simulations in IXu et al.l (120101 l201ll ). The 
AMR is applied only in a volume of (~ 43 Mpc)^ where the galaxy cluster forms near the 
center of the simulation domain. The AMR criteria in this simulation are the same as in 



Xu et al.l (120101 ). During the cluster formation but before the magnetic fields are injected, 
the refinement is controlled by baryon and dark matter density. After magnetic field injec- 
tions, all the regions where magnetic field strengths are higher than 5 x 10"^ G are refined 
to the highest level, in addition to the density refinement. 



The magnetic fields are injected using the s ame method us ed in lXu et al.l ( 120081 . 120091 ) as 
the original magnetic tower model proposed by lLi et al.l (120061 ). We have the magnetic fields 
injected at redshift z = 2.5 in the 30 most massive proto-clusters in the cluster forming region 
of ~ (30 Mpc)^, assuming that there is a super massive blackhole (SMBH) in each of these 
proto-clusters. Out of the 30 proto-clusters with injected magnetic fields, 26 finally form the 
major part of the cluster, while the other four are still in the filaments. We assume that 
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the magnetic fields are from ~ 10^ M© SMBHs witli ~ 1% of outburst energy in magnetic 
form. So we put ~ 1.3 x 10^^ erg magnetic energy into the IGM from each injection, and 
the total injected magnetic energy from all AGNs is about 4 x 10^" erg. The number and 
masses of supermassive black holes used in this simu lation are consistent with black hole 
mass density obtained by Sloan Digital Sky Survey (lYu fc Tremaind |2002| ) . We chose to 
inject the magnetic fiel ds at the redshift z = 2.5, which is close to t he peak of corn oving 
quasar number density (IFan et al.l 1200 ll : iFanI |2006| ) . Previous studies (IXu et al.l 12010) have 
shown that the distribution of the ICM magnetic fields at low redshifts is not very sensitive 
to the exact injection redshifts and injected magnetic energies. 



2.2. Properties of the Galaxy Cluster and Distribution of Magnetic Fields 

Before we present the comparisons between observations and simulations, we first de- 
scribe the properties and the magnetic field distribution of this simulated cluster with multi- 
AGN magnetic field injections. This simulated cluster is a massive cluster with its basic 
properties at redshift z = as follows: Rviriai = 1-66 Mpc, Myiriaiitotal) = 5.68 x 10"*^^ Mq, 
^viriaiig^s) = 7.75 X 10"'^'^ Mq, and Tviriai = 4.52 keV. The magnetic fields in t his cluster 



are sim ilarly distributed as in the clusters of single AGN injection simulations in IXu et al. 



(120 111 ). In Figured], we show the images of the isocontours of the baryon density and mag- 
netic field strength, and a sample of magnetic field lines of this simulated galaxy cluster at 
redshift z = 0. This cluster is disturbed by a big merger after redshift z = 0.25 and is still 
unrelaxed at the current epoch, so it looks irregular in shape. The isocontours of baryon 
density show shell structure visually with the highest baryon densities at the cluster center 
and then the densities decreasing monotonically with radius in all directions. But there is 
no such shell structure of the isocontours of magnetic field strength, showing that magnetic 
fields have no clear dependence on the distance to the cluster center (or the gas density) in 
3-dimensional spatial distribution. This is likely because large amount of magnetic fields are 
generated by the dynamo process of turbulence in additional to the simple compression of 
"frozen-in" magnetic fields. The plot of magnetic field lines further shows that the magnetic 
fields are highly entangled as the magnetic field lines are stretched and twisted by the ICM 
turbulence, while there are few field lines extended out of the cluster. 



The magnetic energy evolution is shown in the top panel of Figure [21 At redshift z = 0, 
the magnetic energy inside the cluster virial radius is 2.37 x 10^° erg, while the total magnetic 
energy in the simulation domain is 3.07 x 10®" erg. So ~ 7 7% of the total m agnetic energy is 



inside the galaxy cluster. The a parameter in Equation 1 of lXu et al.l (l201ll ). which measures 



the relation between magnetic energy and cluster mass, is 1.21, slightly bigger than those 
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1) of the single AGN injection clusters in IXu et al.l ( 1201 ll ). So this cluster, though having 
much more initial magnetic fields, still follows the E^^ oc M^^^.^^ scaling found in single AGN 
model simulations. We also show the spherically averaged radial profiles of the magnetic 
field strengths at z = 0.5, 0.25, and 0.0 in the bottom panel of Figure O Though the 3- 
dimensional dependence of magnetic field strength on radius i s wea k, the averaged magnetic 
field strengths drop with radius monotonically (see IXu et al.l (120111 ) for a detailed discussion 
on the dependence of field strength on radius and gas density). The mean magnetic field 
strengths at z = are just above 1 fiG in the cluster center and decrease slowly with radius 
to about 0.2 fiG at the virial radius. The magnetic field strengths have little evolution 
between z = 0.5 and z = 0, except some bumps related to mergers at z = 0.5. Based on the 
mean magnetic field profiles of clusters in Xu et al. (2011), one would expect an unrelaxed 
cluster such as this to have an irregular radial profile. However, the radial profiles here are 
smooth, most likely due to the fairly regular distribution of the magnetic injection sites. In 
addition, these profiles show that the magnetic field strength does not peak at t he cluster 



center. Since magnetic fields are significantly amplified by the ICM turbulence (jXu et al. 



20101), the density and magnetic field strength peaks do not necessarily coincide. 



3. Comparisons of Simulated Results with Radio Observations 



To see how well the magnetic fields in our simulated clusters fit the observations, we 
compare the synthetic radio observations from our simulations with observations of RM and 
radio halos. In addition to our multi-AGN run, we also include a set of twelve galaxy clusters 
with a wide range of sizes and formation histories with single AGN magnetic field injection 
from lXu et all (|201l[ l. 



3.1. Simulated and Synthetic Faraday Rotation measure Images 

The simulated X-ray and RM images for the multi-AGN run at redshift z = 0.5, 0.25, 
and are shown in Figure |3l The X-ray images are calculated within a box of (4 h~^ Mpc)^ 
(comoving) in the band of 0.1 to 2.4 keV using yt. The simulated RM maps are calculated 
within the same box assuming that a polarized synchrotron source is behinci^ the whole 



^Note that in this way we are overestimating the RM by approximately a factor of V2 with respect to 
the case in which the source is located half-way at the center of the cluster. 
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RMsim = 812 / rieBydy (rad/m^ 



volume 



where rig is the electron number density in cm~'^, By the magnetic field component along 
the line-of-sight in /iG, and L is the size of the computational box in kpc. There is clear 
correlation between the X-ray luminosity and the RM, as the significant RM (> 100 rad 



m ) only appear within bright X-ray (> 1 x 10 erg cm 



-2 s-'] 



regions. The RM patterns 



are patchy and show both large and small scale structures, as expected given the turbulent 
state of the simulated magnetic fields in the ICM. 



In the following, we compare the simulated 



Clarke et al. 


(2001 


) and 


Govoni et al. 


(2010 


)• 



'araday r o tation from the multi-AGN in- 
with observations from 



Xu et al 



(H) 



RM for radio so urces at different irn pact parameters from the centers of a sample of nearby 
galaxy clusters. I Govoni et al.l ( l2010l ) report the dispersion of the Faraday rotation {aRu) cal- 
culated from the distribution of all pixel values of spatially resolved RM images of a sample of 
radio sources in (or behind) nearby galaxy clust ers. Both these indicators derive ultimately 
from EqUl but according to iMurgia et al.l ( 12004 ) the two quantities trace the magnetic field 
on different scales. For an isotropic turbulence, the \RM\ offset from zero would be due to 
ICM magnetic field fiuctuations on scales larger than the projected radio source's size at 
the distance of the cluster. The a^M instead would be due to magnetic field fiuctuations 
on scales smaller than the projected radio source's size. Indeed, the ratio of \RM\/ a^M 
is directly related to the intra-cluster magnetic field power spectrum. In general the \RM\ 
from point-to-point in a cluster is more scattered since it involves a smaller number of mag- 
netic field fiuctuations along the line-of-sight and it is sensitive to the foreground Faraday 
rotation of our Galaxy, which must be subtracted from the data. The a^M is statistically 
more stable, since it depends on a large number of small scale fiuctuations, and it is mostly 
insensitive to unrelated foreground Faraday screens, but it requires high quality RM images 
to be measured. Usually, only a small number of sources are observed per cluster, thus 
only a few limited patches of the magneto-ionic medium are accessible through the observed 
RM images. Indeed, for a proper comparison with observations, the simulated RM images 
have been filtered as typical VLA polarization observations, as described in Appendix A. We 
then calculate the averaged radial profiles of the absolute value of \RM\ and the dispersion 
(jRM- These radial profiles have been obtained from the average of the scatter of the RM 
statistics computed over boxes of 100 x 100 kpc^ which is the typical size of observed RM 
images of cluster radio sources (see Appendix for details). In the following, we refer to the 
filtered RM images as the "synthetic" RM images, to distinguish them from the "simulated" 
RM images shown in the bottom panels of Fig. [31 In the following we perform a qualitative 
comparison of the synthetic RM with the observations with the aim to determine whether 
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the magnetic fields in our simulated clusters are able to provide, at least at the qualitative 
level, a consistent description of the observed magnitude and trends of the RM signal. A 
detailed quantitative analysis would require a much larger statistics of both simulated and 
observed clusters and it is beyond the scope of this work. 



3.1.1. Comparison with Observed \RM\ Profiles 



We first compare the radial profiles of \RM\ of the synthetic RM rn aps from the multi- 
AGN injection simulatio ns, as well as the results from IXu et al.l (1201 ll ). with observations 
from Iciarke et all J200lh . 



We plot the radial profiles of \RM\ of the multi-AGN case at various redshifts in the 
top panel in Figure HI where the lines represent the average of the \RM\ of all the boxes at 
a given radius from the cluster center. It must be noted that, due to the fluctuations of the 
intra-cluster magnetic field on scales larger than 100 kpc (i.e. the size of the boxes in which 
the RM statistics is calculated), the synthetic profiles are characterized by a strong intrinsic 
scatter, with a dispersion of almost one order of magnitude around the average. We show 
the rms scatter as a shaded region for the RM profile at redshift 2; = 0. The \RM\ profiles 
at z = 0.25 and z = 0.5 have a similar scatter (not shown for graphical simplification). 

Then we show the profiles from the relaxed and unrelaxed cluster^ in Xu et al. (2011) 
at 2 = in the middle and botto m panels, respectively. These clusters are called by the 



same names as m 



Xu et al.l (I2OIII ). There are t wo profiles of clu ster Rl shown, designated 



as Rla and Rib, as the simulations A and B in IXu et al.l (120101 ). which are from the same 
simulated cluster with two different amounts of injected magnetic fields. All the simulated 
clusters will be named in this way throughout the paper. The temperatures of these clusters 
are 7.65, 5.9, 4.14, 3.04, 2.16, and 1.41 keV, for the relaxed clusters from Rl to R6, and 
10.26, 4.84, 4.78, 4.59, 4.34, and 3.43 keV, for the unrelaxed clusters from Ul to U6. The 
shaded regions in middle and bottom panels correspond to the intrinsic scatter of the Rla 
and Ula simulations, respectively. 



The results of observations from lClarke et al.l (120011 ) are over plotted. The data are from 
different clusters, since it is still not possible to measure a large number of RM at different 
radii in a single cluster due to the limited sensitivity of current polarization observations. 
Two different symbols are used to show whether the cluster temperatures are hotter (filled 



■^Note that we name "relaxed" those simulated clusters which are in the post merger stage. These systems 
should not be confused with the relaxed cool-core clusters often referred in the literature. 
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dots) or cooler ( empty dots) than 5 keV. The temperatures of the observe d clusters are taken 
from literature JProust et aliboosl : Ichen et al.lboOTt ISogdan et aPboill ). 



The observed \RM\ from lClarke et al.l (120011 ) have a clear trend that the absolute values 
of RM drop gradually with increasing radius, while RM have typical values of 100s rad m~^ 
near the cluster centers. The synthetic \RM\ profiles from all simulations clearly show similar 
trends, but have some different features, likely due to the differences in cluster properties 
and in the cluster formation histories. The \RM\ profiles from the multi-AGN simulation 
are smooth and match the small \RM\ from observations. It is likely because this simulated 
cluster is not very large, just a little cooler than 5 keV, and its magnetic volume filling is 
high. Furthermore, the \RM\ profiles, similar to the magnetic field strength profiles, have 
little evolution between z = 0.5 and z = 0. The profiles from single AGN injected clusters 
clearly show a positive trend between the cluster mass (temperature) and the |-RM|, and 
the profiles from large clusters cover those large observed |i?M|. The \RM\ from relaxed 
clusters are usually larger than the \RM\ from unrelaxed clusters in similar size. 

We note once again that the \RM\ is an intrinsically scattered indicator. The nature of 
the scatter is not due to instrumental errors but rather represents an intrinsic characteristic 
of the intra-cluster magnetic fields, as illustrated by the shaded regions shown in Figure [Hand 
by the spread of observed data points as well. The scatter of the \RM\ inside a single cluster 
is comparable to the difference between the average \RM\ profiles among different clusters. 
It is then obvious that a detailed comparison would required a much larger statistics of both 
simulated and observed galaxy clusters. Nevertheless, from this first comparison, we can 
conclude that most of the simulated \RM\ profiles lie within the range of values spanned by 
the data, with an overall agreement which is within the expected intrinsic scatter. There are, 
however, a few exceptions like simulations R6 and U2b whose \RM\ profiles are well below 
the data. On other hand, this is not surprising since th ese are arnong th e least magnetized 
systems in the sample of simulated clusters described in IXu et al.l (120111 ). 



3.1.2. Comparison with Observed a hm Profiles 



We plot the radial profiles of crm of the multi-AGN simulation at redshifts z = 0.5, 
0.25, and in the top panel of Fig ured while the profiles of relaxed clusters and unrelaxed 
clusters at z = in lXu et al.l (120111) are shown in th e middle and bottom panels, respectively. 
The results from observations (iGovoni et al.ll2010l ) are over plotted, ar id are represe n ted b y 
different symbols in the cluster temperature, as listed in the table 5 of iGovoni et al. J2010h . 
As in Figure 4, we show, as shaded regions, the intrinsic scatter of the a^M profiles for the 
multi-AGN simulation at redshift z = and for the Rla and Ula simulations. 
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The apiM have similar radial profiles as |-RM|, but data and synthetic profiles (as illus- 
trated by the shaded regions) are much less dispersed. The (Trm near the cluster centers 
are of the order of a few 100 rad/ut^ and gradually drop with radius. The flattening seen 
at large radii is an instrumental effect: the minimum a^M that can be measured is limited 
by the noise in the observed RM images and we account for this effect in the synthetic RM 
images (see Appendix A). 

The synthetic a^M profiles are generally within the range spanned by the data. The 
profiles for unrelaxed clusters are bumpy and the scatter is somewhat larger, probably 
as a result of recent mergers occurred in these systems. In addition, the synthetic profiles 
reproduce the statistical treri d observed in the data between the grm and the cluster tem- 
perature (iGovoni et al.ll2010l ). Namely, the more massive, hotter galaxy clusters are also 
characterized by higher values of auM- Indee d, simulations expect that hotter clusters to 
be also more magnetized (see IXu et al.ll201l[ ) and, as a co n seque nce, they should present 
also higher aRM- However, as pointed out by iGovoni et al.l (l2010l ). hotter clusters are also 
generally denser, at least in their sample, and thus it is not easy to determine which is the 
dominant effect between magnetic field strength and gas density. This holds in general for 
all simulated aRu profiles. In fact, aRM basically scales as a function of gas density and 
magnetic field strength as 

-, 0.5 



aRM oc 812 



A. 



dy 



(2) 



wher e A^ is the auto-correlation length of the magnetic field fluctuations (see e.g. lMurgia et al, 
20041 ). The value of aRM depends directly on the product of B and rig. From Figure Owe 
show that simulations reproduce quite well the observed radial profiles of aRM, but we are 
not able to tell if this is due to the right combination of B and Ug. 



3.1.3. Comparison with Observed aRM - Sx Relation 



A possibility to break the degeneracy between magnetic field strength and gas density 
is to investigate the rel ation ship between ctrm and the X-ray surface brightness Sx- In- 
deed, lDolag_et_alJ (I2OOII ) and lGovoni et al.l (I2OIOI ) showed that there is a positive correlation 
between these quantities. In fact, 



Sx oc 



n. 



'Vrdy. 



(3) 



Thus, by comparing aRM and Sx we can try to reduce the degeneracy between and B in 
Eq.El 
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We present the (Trm ^ Sx relations from our simulations and the observations from 



Govoni et al.l ( 120101 ) in Figure |6l The results from the multi-AGN simulation at redshifts 



z = 0.5, 0.25, and are plotte d in the top panel, while the results from relaxed and unrelaxed 



clusters from IXu et al.l ( 120111 ) at redshift z = are shown in the middle and bottom panels, 
respectively. The aRM and the corresponding X-ray surface brightnesses have been computed 
consistently from the same regions of 100 x 100 kpc^. No redshift corrections of a^M and 
Sx had been done, since all the observations are at low redshifts (< 0.1) and the corrections 
are not significant. As in Figure 4 and 5, the shaded regions represent the intrinsic scatter 
of aRM for the multi-AGN, Rla, and Ula simulated clusters. 

We note that the simulated profiles show the same trend as data, in particular the 
slope of the simulated a^M ^ relation is close to the observed one with a clear power-law 
regime at high a^M and Sx- At very low a^M and Sx the synthetic correlations flatten 
because the cluster RM signal fall below the detection threshold of the RM observations, 
which is about 30 rad/m^. We recall that this is one of the most relevant features introduced 
by the instrumental filtering (see Appendix A). For the unfiltered simulated clusters (not 
shown), the power law regime of the ctrm ^ Sx relation extends for more than three order 
of magnitudes down to very low X-ray surface brightness and RM. 

Although the synthetic profiles show a slope consistent with the observed trend, they 
appear to be systematically shifted toward high Sx- This suggests that the simulated clusters 
could be slightly over-dense or under magnetized than the observed ones. Since Sx also 
increases with the square root of T, it could be also possible that simulated clusters are 
systematically hotter than the observed ones. However, we checked that the distributions 
of global temperature for simulated and observed clusters are quite similar. Indeed, the 
average temperature for simulated cluster is of 4.7 keV which is even lower than the average 
temperature of 6.0 keV found for the galaxy clusters in Govoni et al (2010). This shift is 
more evident for the relaxed and unrelaxed simulated clusters with a single AGN injection, 
as the under magnetization is likely. For example, the radial profile of simulated cluster Rla 
reproduces the upper bound of data in Figure [5] but lies just in the middle of the ctrm - Sx 
relation and the simulated Sx peak exceeds that of data by roughly a factor of four. 

In summary, the comparison of Figures O and M suggests t hat simulated cluste rs should 



by slightly less dense and more magnetized to fit those in the iGovoni et al.l (120101 ) sample. 
If simulated magnetic fields were a factor of two stronger while gas densities a factor of two 
lower, ajiM would remain unchanged but the expected a^M ^ Sx relations would shift to 
the left by a factor of four. The multi-AGN injection simulation seems to simultaneously 
provide a consistent description to both the aRM radial profiles and the aRM - Sx relation, 
suggesting a better combination of B and Ue for this cluster. Again, there is no evolution of 
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these relations from z = 0.5 to 2; = for the multi-AGN injection modeL It maybe because 
this cluster is well magnetized during all this time. 



3.2. Simulated and Synthetic Radio Halo Images 



The comparison between observed and synthe tic radio halo imag e s is another tool to 



investigate magne tic fields in clusters of galaxies (jMurgia et al. 2004 : Govoni et al 



2006 



Vacca et al.ll2010l ). This approach consists of simulations of radio halo images obtained by 
illuminating cluster magnetic fields models with a population of relativistic electrons. 

In the specific case considered in this work, the radio halo modeling simulations have 
been performed at 1.4 GHz with a bandwidth of 25 MHz. At each point of the compu- 
tation grid, we calculate the synchrotron emissivity by convolving the emission spectrum 
of a single relativistic electron with the particle energy distribution of an isotropic popula- 
tion of relativistic electrons whose distribution follows N{'-y,9) = KQ'-y~^{sm9)/2, where 7 
is the electron's Lorentz factor, while 6 is the pitch angle between the electron's velocity 
and the local direction of the magnetic field. The energy density of the relativistic electrons 
is Uei = rrieC^ pmi^'^ NM'jd'j, where 7min and 7niax are the low and high energy cut offs of 

Tmin 

the energy spectrum, respectively. In our modeling we assume equipartition between the 
magnetic field energy density mb = B^/Stt and u^i at every point in the cluster, therefore 
both energy densities have the same radial decrease. Radio halos have a typical spectral 
index a = 1.3, thus we adopt an electron energy spectral index 6 = 2a + 1 = 3.6. The 
low-energy cutoff ■y^nin and normalization Kq are adjusted to guarantee u^i = ub, in practice 
we fixed 7min=300 and let Kq vary. The energy losses o f the relativist ic electrons in the ICM 
reach a minimum for this value of the Lorentz factor ( ISaraziru 119991 ). For the steep power 
law energy spectrum we considered, the exact value of 'jmax has only a marginal impact on 
the modeling, we adopted a fixed value of •jmax = 1.5 x 10"^ for consist ency with the radia l 
steepening observed in the spectral index images of some radio halos (IFeretti et al.ll2012l ). 
Most observations have failed to detect polarization in radio halos, and thus in the present 
paper we focus on the analysis of the total intensity emission. Nevertheless, the investigation 
of the polarized intensity of the simulated radio halos is important and will be the subject 
of a next publication. 

We produced radio halo images for the simulated cluster (at redsh ift z=0) with multi- 
AGN injections, as well as for the simulated clusters by lXu et al.l (l201l[ ): Rla, Rib, R2, R3, 
R4, R5, R6, Ula, Ulb, U2a, U2b, U3, U4, U5, U6. For a proper comparison with observa- 
tions, the images have been filtered as typical VLA radio halo observations, as described in 
Appendix B. For each cluster we obtained a synthetic radio image with an angular resolution 
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of 50" FWHM and a noise level of 0.1 mJy/beam (la). We found that in multi-AGN, Rla, 
Rib, R2, R3, Ula, Ulb, and U4 the radio signal is above 3cr and the radio emission is strong 
enough to be detected in a typical VLA observation. In all the other simulations the signal 
is below the 3cr limit and only an upper limit to the flux density can be given. Following 
the same approach adopted for the RM simulations, we compared the synthetic radio halos 
with the observations. In particular, we focus on the offset between the radio and the X-ray 
peaks, on the correlation between the global radio luminosity at 1.4 GHz and the X-ray 
luminosity (Pi.4ghz - Lx), and on the correlation between the Largest Linear Size (LLS), 
as measured from the 3cr isophote on the synthetic images, and the radio power. These 
properties are deriv ed from the synthetic halo images and are compared to the data taken 
from the literature f lFeretti et al.ll2012t iGovoni et al.ll2012[ ). The evolution of the relativistic 
electron component is not treated explicitly in our simulations and hence we cannot predict 
the formation of radio halos. Our aim here is to check whether, given a reasonable energy 
spectrum for the synchrotron electrons, the simulated magnetic field produces radio halos 
whose global properties are in line with the observations. 



3.2.1. Comparison with the Observed Radio to X-ray Offsets 



First we compare the offset between the X-ray and the radio peaks (see iFeretti et al. 



201ll ). This is a quite easy observable to analyze and it can be related to the dynamic 
state of the ICM. In particular, it could be expected that unrelaxed clusters with ongoing 
mergers will show larger offsets with respect to relaxed clusters. Mor eover, this offset co uld be 
directly related to the magne tic field power spectrum on large scales (IVacca et al.ll2010l ). Very 
recently, iGovoni et al.l (120121 ) measured this offset for a sample of galaxy clusters containing 
radio halos by analyzing VLA and ROSAT data. In Figure [7] we show the radio contours 
of the detected synthetic halos overlaid on the simulated cluster X-ray emission in the 0.1- 
2.4 keV band. We smoothed the X-ray emission to a resolution of 45" in order to mimic 
a ROSAT PSPC observation. In this way we can compute the offset between the X-ray 
and the radio peaks in a similar way as in most data of galaxy clusters hosting radio halos. 
We note that the synthetic radio halos can be quite asymmetric with respect to the X-ray 
gas distribution. In Figu re [8] we compare th e offset between the radio and X-ray peaks for 
a sample of radio halos ( iGovoni et al.l 120121 ) with simulations. The offset measured in the 
synthetic radio images are on the range of the observed ones, although a detailed comparison 
of the two distributions would require a much larger statistics for both simulations and data. 
We note that the offset is comparatively small for the multi-AGN injection and for the 
relaxed clusters Rla, Rib, R2 and R3. On the other hand, the offset is particularly large for 
the unrelaxed clusters Ula and Ulb, as expected. In fact, the offset is so extreme for these 
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system s that they could be even classified as cluster roundish relics according to lFeretti et al. 
fl2012h . 



3.2.2. Comparison with Observed Pi ^q^^^Lx and Piaguz^LLS Relations 

In Figure M we plot the radio power calculated at 1.4 GHz versus the cluster X-ray 
luminosity in the 0.1—2.4 keV band and versu s the Largest Linear Size. Blue dots r efer to 
the data published in the literature (taken from lFeretti et al.ll2012l : lGovoni et al.ll2012[ ). while 
red triangles refer to the simulations. In the simulations Lx has been calculated within a 
circle of 1 Mpc in radius, while Pi.4ghz and LLS have been obtained by considering the radio 
halo emission up to the 3a level. Eight out of 16 simulated clusters present a detectable radio 
halo under the hypothesis of equipartition condition. Among the "radio loud" clusters we 
have the relaxed clusters Rla and Rib and the unrelaxed clusters Ula and Ulb. These four 
clusters are also the most X-ray luminous systems of the simulation set and they all have 
about the same X-ray luminosity of ~ 10"^^ erg/s, but are quite different in radio luminosity. 
In fact, the diffuse radio emission in Rla and Rib is about one order of magnitude more 
powerful than in Ula and Ulb. The differen ce could b e relate d to their different evolutionary 
histories. According to the classification in IXu et al.l ( 1201 ll ). R-type clusters have accreted 
half of their actual masses hj z = 0.5 while U-type clusters gained more than half of their 
final masses from z = 0.5 to 2; = 0. The decaying turbulence originated from the early 
mergers in Rla and Rib has had enough time to diffuse and amplify the magnetic field over 
a large portion of the cluster volume, which results in a powerful and extended radio halo at 
the cluster center. On the other hand, the merging process is still in progress for Ula and 
Ulb, and consequently there has not been enough time to amplify the magnetic field, to the 
level found in Rla and Rib. The synthetic radio halos in Rla, Rib, Ula, Ulb, together with 
those in R2, R3, U4 agree with the observed correlations between the radio power versus 
the cluster X-ray luminosity and between the radio power versus the radio halo size. The 
synthetic radio halo in the case of multi-AGN magnetic field injections lies perfectly in the 
relation between radio power versus the radio halo size, but is over-luminous in radio with 
respect to the clusters R3 and U4 which have a very similar X-ray luminosity. This may 
imply that clusters in which the injection of seed magnetic field by AGNs has been particular 
efficient may also host radio halos with a power higher than average. We speculate that this 
is could be a possible expla nation also for the few radio ove r-luminous outliers kno wn so 
far in the literature: A523 (IGiovannini et al.l I2OIII ). A1213 (IGiovannini et al.ll2009l ). and 
0217-1-70 (IBrown et al.l 1201 ll ). For the remaining 8 simulated clusters whose radio emission 
is below the detection threshold we can only derive upper limits. In order to calculate the 
upper limit on total radio power Pi.4ghz for the undetected radio halos we must first obtain 
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an estimate of their putative size. On the basis of the Lx-LLS correlation (not shown), 
we extrapolated the putative radio halo size from its X-ray luminosity. The upper-limits on 
total radio power Pi.4ghz were then calculated with the assumption that the average surface 
brightness of the diffuse emission over the putative halo size is lower than the Sex noise level 
of the corresponding radio image. We note that the clusters with a non-detected radio halo 
emission are the fainter X-ray clusters. 

In summary, the magnetic field models presented here seem to generate halos of the right 
power and size, under the hypothesis of energy equipartition. By reversing the argument, 
we can deduce that the central magnetic field strength and radial decline in our simulations 
combine to produce a volume average magnetic field that matches the equipartition magnetic 
field estimates for observed radio halos. 



Conclusions 



To demonstrate the effects of location and number of injections on the ICM magnetic 
fields, we study the magnetic field evolution of a galaxy cluster with magnetic fields injected 
by as many as 30 AGNs at redshift z = 2.5. The cluster is well magnetized at low redshifts, 
the final mean cluster magnetic fields are ~ 1 fiG at the cluster center and ~ 0.2 fiG at its 
virial radius. Though having wider spreading and higher volume filling, the magnetic field 
distribution of this cluster is not significantly different from those of our previously single 
AGN injection simulations. The additional injected magnetic fields do not increase the final 
magnetic energy and field strength dramatically, but only increase the volume filling of the 
cluster, especially in the outer parts of the cluster. The total magnetic energy of this cluster 
is consistent with the scaling relation with the square of cluster mass in the high end. We 
also find smoother radial profiles of magnetic field strength after major mergers, likely due 
to the fast amplification of the high volume filled magnetic fields. 

We produce VLA filtered synthetic Faraday rotation images from a set of simulated 
galaxy clusters and compare them with real observations. The radial pro files of the synthetic 
\R M\ and aRM from sini ulations agree with the observation results from lClarke et al.l (120011 ) 
and iGovoni et al.l ( l2010l ). The \RM\ and aRM, just like magnetic field strengths, are high 
in large (hot) clusters, as suggested by observations ( iGovoni et al.ll2010l ). We also compare 
the ajiM - Sx relation between synthetic results and observations, and show that they are 
in a broad agreement, although there is an indication that the simulated clusters could be 
slightl y over-dense and les s magnetized (up to about a factor of two) with respect to those 
in the 



Govoni et al.l (|2010[ ) sample. 
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We generate synthetic radio halos by assuming the energy equipartition between mag- 
netic fields and non-thermal electrons, and we compare the global properties of the mock 
radio halos with the real ones. The offset between the radio and the X-ray peaks for the syn- 
thetic halos are in the range of the observed values recently measured in a sample of clusters 
hosting radio halos. In particular, we find the offset is small in relaxed clusters and large in 
unrelaxed clusters, suggesting that this indicator is strongly related to the dynamic state of 
these systems. Furthermore, the synthetic radio halos agree with the observed correlations 
between the radio power versus the cluster X-ray luminosity and between the radio power 
versus the radio halo size. 

Therefore, we can draw the important conclusion that the combined effects of the central 
magnetic field strength and the radial decline in our simulations can produce a volume 
average magnetic field that matches the equipartition magnetic field estimates for observed 
radio halos and, at the same time, is able to explain the observed valu e of Faraday rotation 



me asure of radio sources in clusters. This results confirm the findings by lMurgia et al.l (120041 ) 



and lGovoni et al.l (120061 ) that if the ICM magnetic fields fiuctuate over a wide range of spatial 
scales and decline in strength with radius, it is possible to reconcile the equipartition magnetic 
fields in radio halos and the magnetic field derived from the RM in galaxy clusters. 
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A. RM Images Filtering 

When comparing the expectation of numerical simulations of cluster Faraday rotation 
images with data one must take into account of the instrumental filtering affecting the 
observed RM images. The observed RM images of cluster radio galaxies considered in this 
work have been obtained through the well known A^— law relating the observed polarization 
position angle at wavelength X = c/u to the intrinsic polarization angle \&o ^ given 
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position in the radio source: 



^^/{\) = ^o + RMxy{l + zy. (Ai) 



By observing the radio source at different frequencies it is possible to obtain both RM 
and \l/o by a hnear fit to Eq. lAll The approach is straightforward in principle but compli- 
cated in practice by a number of observational issues. One is the beam depolarization. The 
spatial distribution of polarization angles is not uniform across the radio source both because 
of the intrinsic morphology of the source itself and because of the Faraday rotation. Indeed, 
the incoherent sum of the disordered polarization vectors inside the observing beam results 
in an attenuation of the polarized signal. Another effect to consider is the bandwidth de- 
polarization. The Faraday rotation of the polarization angle inside the observing frequency 
bandwidth leads again to an attenuation of the signal. Finally, another difficulty to face 
is related to the n-vr ambiguity of the observed polarization angles which must be resolved 
by the fitting algorithm. All these effects are more prominent at low frequencies and in 
general lead to non-linear distortions of the A^— law which are fairly a bit more complicated 
than the deviations induced by a simple propagation of the error measurements from the 
polarized intensity images. The consequent error on the observed RM is difficult to quantify 
analytically. 

For all these reasons, we filtered the simulated Faraday rotation images before comparing 
them to data. The filtering is performed with the FARADAY code and consists of the 
following steps: 

i) first, we created two mock images of 4096 x 4096 pixel in size for the distribution of 
intrinsic polarization intensity and angle at a resolution eight times finer than that of 
the simulated Faraday rotation images (i.e. ^ 1.3 kpc). To mimic the intrinsic vari- 
ations seen in real sources, the polarized intensity and angle distributions are created 
randomly with Gaussian fluctuations on spatial scales smaller than 100 kpc. 

ii) the simulated RM images are gridded to the size of the intrinsic polarization images and 
we employed Eq. lAll to rotate the intrinsic polarization vectors at freq uencies of 8465, 



8085, 4935, 4735, and 4535 MHz. This is the frequency setup used by iBonafede et al. 



teOlOt l for the RM study of the Coma cluster, but it can be considered representative 



also of most other similar studies based on VLA observations; 

iii) we produced synthetic polarization images of Stokes parameter (/{u) and Qii^) at each 
frequency by adding a Gaussian noise such that we have an average signal-to-noise ratio 
of (P/ap) = 5 for the polarization intensity a 8465 MHz. The polarized intensity at 
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lower frequencies is scaled with a spectral index of 0.8, i.e. -P(z^) = P8465('^/8465) 
so that the signal-to-noise ratio is slightly higher at lower frequencies; 

iv) to simulate the beam depolarization the synthetic (/{u) and Qi^i') images are convolved 
with a FWHM beam of 4 x 4 pixels (which corresponds to about 3x3 arcsec at z = 0.1). 
Moreover, to simulate bandwidth depolarization, the synthetic polarization images are 
averaged over a bandwidth of 50 MHz centered on each frequency; 

v) finally, the synthetic and Q(z^) images have been gridded to the original size 
and resolution of the simulated Faraday rotation image and then analyzed as if they 
were real polarization observations. The synthetic RM images are created pixel by 
pixel by fitting the "observed" polarization angle versus the squared wavelength for all 
the frequencies. To reduce the problems associated with n-vr ambiguities, the fitting 
algorithm in FARADAY can perform a sequence of improvement iterations. In the first 
iteration, only a subset of high signal-to-noise pixels is considered. In the successive 
iterations, lower signal-to-noise pixels are gradually included and the information from 
the previous iteration is used to assist the fit of the A^-law. In the procedure we clipped 
all pixes where error on the polarization angle at any frequency was greater than 15° 
and those where the reduced of the fit was worse than 10, as usually done with real 
data. 

A comparison of simulated and synthetic RM images is shown in the top panels of 
Figure [10] for the case of the multi-AGN injection simulation at 2; = 0. The instrumental 
filtering results in a RM noise of about 30 rad/m^. In the central region of the cluster, 
where the RM signal is strong, the two images are very similar, although the fit failed in a 
few isolated pixels (gray pixels). However, in the peripheral regions of the cluster the RM 
signal is too weak to be detected with this observational setup which is based on relatively 
high-frequencies. Middle panels of Figure [10] show the fit of the A^-law (solid lines) at the 
two different locations indicated in the synthetic RM image. The dotted lines correspond 
to the trend expected without instrumental filtering. Although the simulated and synthetic 
RM can be different from pixel to pixel, their statistical properties are quite stable. In 
the bottom panels of Figure [10] we show a comparison of the simulated and synthetic RM 
profiles of \RM\ and a^M as a function of the distance from the cluster center. The RM 
statistics is calculated over boxes of 100 x 100 kpc^ in size, to reproduce the typical size 
of observed RM images of cluster radio sources. The profiles shown in FigHD] represent a 
radial average obtained by exponentially smoothing with a length-scale of 50 kpc the \RM\ 
and aRM calculated in all the boxes. We note that the simulated and synthetic profiles of 
\RM\ match each other, this means that the instrumental filtering has a minor impact on 
this indicator. While the flattening of the synthetic a^M profile at large radii is the most 
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relevant effect introduced by the instrumental filtering. The shaded regions represent the rms 
scatter around the average profiles. This scatter is not originated by instrumental errors but 
rather represents an intrinsic characteristic of the intra-cluster magnetic fields. In particular 
we note that the \RM\ is intrinsically much more scattered than a^M- This is due to the 
fluctuations of the intra-cluster magnetic field on scales larger than 100 kpc, i.e. the size of 
the boxes in which the RM statistics is calculated. 

In Figure [TT] we show the comparison of simulated and the synthetic RM images for the 
case of the simulated cluster Rla. The stronger Faraday rotation gradients in this cluster 
causes large depolarization regions in the center. Close to these regions the fit of the A^-law 
can be troublesome, as shown in the middle-left panel. Here, we show the results for a 
pixel located at the interface of two RM cells of opposite signs. The beam and bandwidth 
depolarization caused by the strong RM gradient results in a deviation from the original A^- 
law (dotted red line). The fitting algorithm then fails to solve for the correct n-vr ambiguities 
and the resulting RM and are completely wrong (continuous blue line and dots). 

However, the radial profiles of the simulated and synthetic \RM\ and a^M shown in the 
bottom panels are in a remarkable agreement. 

B. Radio Halo Images Filtering 

As a consequence of their low surface brightness, radio halos are generally studied 
through interferometric images at a relatively low spatial resolution. Thus, although ra- 
dio halos may have an intrinsic filamentary structure, the details of their morphology are 
often very poorly seen when observed with typical resolutions of 40 — 60". In addition, inter- 
ferometers filter out structures larger than the angular size corresponding to their shortest 
spacing, thus a loss of flux density is expected for radio halos with an extended angular 
size. Indeed, in order to compare the simulated radio halos to the observations we produced 
synthetic radio halo images which include the effects of the instrumental filtering. 

In particular, we performed the following steps: 

i) first we performed mock radio halo images at 1.4 GHz with a bandwidth of 25 MHz 
by integrating the synchrotron emissivity at each point of the computation grid along 
the line-of-sight. 

ii) these full resolution images are mapped as they would appear on the sky at a redshift 
z=0.2, which is the average distance of the known radio halos. We also k-corrected 
their surface brightness for the cosmological dimming; 
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iii) we filtered the resulting images in the Astronomical Image Processing Syste with 
the set-up which has been typically used in many of the pointed interferometric obser- 
vations of radio halos reported in literature. In particular, we converted the images in 
the expected visibility data for a 2 hours time-on-source observation with the VLA in 
D configuration; 

iv) the synthetic visibility data-set was then imaged in AIPS as a real observation. At the 
end, we obtained a filtered synthetic image of the radio halo characterized not only by 
the proper noise and spatial resolution, but also including the subtle effects induced 
by the missing short-spacing and by the imaging algorithm. 

Examples of the radio halo image filtering are presented in Figure [12] for the simulated 
cluster with multi-AGN magnetic field injections at 2; = and for the simulated cluster Rla. 

The full resolution radio halo images are shown on the left column panels of Figure [121 
In these images we can appreciate the full extension of the radio halo emission and also the 
fine details of its intrinsic filamentary structure. 

The synthetic radio halo images after the filtering process are shown on the right column 
panels of Figure [121 The synthetic radio images have an angular resolution of 50" FWHM 
and have a noise level of 0.1 mJy/beam (Icr). Most of the fine details of the halo structure 
are not visible anymore and also a significant fraction of the radio halo results well below the 
detection threshold, especially for the intrinsically fainter radio halo of multi-AGN magnetic 
field injections simulation. 
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Fig. 1. — Isocontours of baryon density (left), magnetic field strength (middle), and magnetic 
field lines at 2; = 0. The field of view is 6 Mpc. The color ranges are 1 x 10~^^ (dark) to 3 
X 10~^^ (bright) g cm~^ for density, and 1 x 10~^° (blue) to 2 x 10~^ (red) G for magnetic 
fields. Different magnetic field lines are represented by different colors. 
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Fig. 2. — Top: Evolution of the total magnetic energy (solid line) in the simulation. The 
magnetic energies inside a box of (4 h~^ Mpc)^ comoving centered at the cluster center and 
inside the virial radius of the cluster are marked at redshifts z = 0.5, 0.25, and 0. At z = 0.5, 
most of the magnetic fields (80% of the total magnetic energy ) still reside in the incoming 
halos, and are out of the major cluster. Bottom: The spherically averaged radial profiles of 
RMS magnetic field strength at redshifts z = 0.5, 0.25, and 0. 
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Fig. 3. — Top panel: Simulated X-ray emissions observed from the y direction in the 0.1 to 
2.4 keV band of the cluster at redshifts z = 0.5, 0.25, and 0. Each image covers an area 
of 4 h-i Mpc X 4 h-i Mpc comoving, while the X-rays are calculated from boxes of (4 h ^ 
Mpc)^. The color range is from 1 x 10~^ (blue) to 4 x 10^^ (red) erg s~^ cm~^. Bottom 
panel: Simulated Faraday rotation measure of the cluster by integrating through the cluster 
on the y direction at z = 0.5, 0.25, and 0. They are calculated from the same boxes as for 
the X-ray. The color range is from —400 (blue) to 400 (red) rad m~^. 
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Fig. 4. — Azimuthally averaged radial profiles of |-RM|. The panels, from top to bottom, 
show the results of the multi-AGN run at different redshifts, a set of relaxed clust ers, and a 
set of unrelaxed clusters, respectively. The observations from I Clarke et al.l (120011 ) are over- 
plotted in two groups based on whether the clusters are hotter (filled dots) or cooler (empty 
dots) than 5 keV. The shaded regions represent the rms scatter around the average profiles 
for multi-AGN, Rla, and Ula runs, respectively. 
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Fig. 5. — Radial profiles of the dispersion of RM distribution. The panels, from top to 
bottom, show the results of the multi-AGN run at different redshifts, a set of relaxed clusters. 



and a set of unrelaxed clusters, respectively. The observational data from iGovoni et al. 



(I2OIOI ) are over-plotted in two groups based on whether the clusters are hotter (filled dots) 
or cooler (empty dots) than 5 keV. The shaded regions represent the rms scatter around the 
average profiles for multi-AGN, Rla, and Ula runs, respectively. 
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Fig. 6. — Dispersion of RM distribution as a function of the X-ray surface brightness. The 
panels, from top to bottom, show the simulated results of the multi-AGN run at different 
redshifts, a set of rel axed clusters, and a set of unrelaxed clusters, respectively. The ob- 
servational data from iGovoni et al.l (120 lOf ) are over-plotted in two groups based on whether 
the clusters are hotter (filled dots) or cooler (empty dots) than 5 keV. The shaded regions 
represent the rms scatter around the average profiles for multi-AGN, Rla, and Ula runs, 
respectively. 
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Fig. 7.— Simulated clusters: multi-AGN, Rla, Rib, R2, R3, Ula, Ulb, U3, and U4. Radio 
halo emissions mapped to sky and filtered like typical VLA observations in D configuration 
with an integration time of 2 hours per cluster. The radio images have an angular resolu- 
tion of 50" FWHM. The contour levels start at 0.3 mJy/beam {3a) and increase by \/2. 
Radio halo contour levels are overlaid to the cluster X-ray emission in the 0.1-2.4 keV band 
smoothed at 45" resolution to mimic a ROSAT PSPC observation. The radio emission from 
U3 is barely detectable, but we consider it as undetected since in a real image it would be 
indistinguishable from a background radio source. 



- 33 - 



10 



Si 

CP 

6 
3 




' I ' I ' I ' I 

I I OBSERVATIONS 
^ SIMULATIONS 



100 200 300 400 500 600 700 800 900 

Radio— Xray Offset [kpc] 



Fig. 8. — Offset between the radio and X-ray peaks for a sample of radio lialos (from Govoni 
et al. 2012), compared witli simulations. 
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Fig. 9. — Left panel: radio power of radio halos at 1.4 GHz versus the cluster X-ray 
luminosity in the 0.1—2.4 keV band. Right panel: radio power of halos at 1.4 GHz versus 
their largest linear size (LLS) measured at the same frequency. Blue dots are observed 
clusters, red triangles are simulated clusters, while arrows indicate upper limits on the radio 
power of simulated clusters. The LLS are measured from the 3cr isophote on the halo images. 
The data are taken from the recent compilations by Feretti et al. (2012) and Govoni et al. 
(2012). 
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Fig. 10. — Simulated and synthetic RM images of multi-AGN injection cluster at z = are 
shown in the top left and right panels, respectively. Middle panels show the fit of the A^-law 
(solid lines) at the two locations indicated in the synthetic RM image. The dotted lines 
correspond to the trend expected without instrumental filtering. The value of the unfiltered 
Faraday rotation value {RMsim) is also reported in each plot. Bottom left and right panel 
show the simulated and synthetic \RM\ and a^M profiles, respectively. The profiles represent 
an exponentially smoothed radial average of the RM statistics calculated in boxes of 100 x 100 
kpc^. The shaded regions represent the rms scatter around the average profiles. 
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Fig. 11. — Same as Figure flOl but for the simulated cluster Rla. The stronger Faraday 
rotation gradients in this cluster causes large depolarization regions in the center (gray 
pixels). Close to these regions the fit of the A^-law can be troublesome, as shown in the 
middle- left panel. However, the radial profiles of the simulated and synthetic \RM\ and ajiM 
shown in the bottom panels are in a remarkable agreement. The shaded regions represent 
the rms scatter around the average profiles. 
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Fig. 12. — Full resolution total intensity radio image at 1.4 GHz, obtained from the simulated 
cluster with multi-AGN magnetic field injections at redshift z = (top-left panel) and the 
simulated cluster Rla (bottom-left panel). The right panels show the synthetic radio halo 
emission mapped to sky and filtered like a typical VLA observation in D configuration with 
an integration time of 2 hours. The synthetic images have a resolution of 50" FWHM beam. 
The contour levels start at 0.3 mJy/beam (3a) and increase by \/2. 



